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Superfluid turbulence is governed by two dimensionless parameters. One 
of them is the intrinsic parameter q which characterizes the relative value 
of the friction force acting on a vortex with respect to the non-dissipative 
forces. The inverse parameter q~ l plays the same role as the Reynolds num- 
ber Re = UR/u in classical hydrodynamics. It marks the transition between 
the "laminar" and turbulent regimes of vortex dynamics. The developed tur- 
bulence, described by a Kolmogorov cascade, occurs when Re 3> 1 in classical 
hydrodynamics. In superfluids, the developed turbulence occurs at q <C 1. 
Another parameter of superfluid turbulence is the superfluid Reynolds num- 
ber Re s = UR/k, which contains the circulation quantum k characterizing 
quantized vorticity in superfluids. The two parameters q and Re s control 
the crossover or transition between two classes of superfluid turbulence: (i) 
the classical regime, where the Kolmogorov cascade (probably modified by 
the non-canonical dissipation due to mutual friction) is effective, vortices 
are locally polarized, and the quantization of vorticity is not important; and 
(ii) the Vinen quantum turbulence where the properties are determined by 
the quantization of vorticity. The phase diagram of these dynamical vortex 
states is suggested. 

PACS numbers: 43.37.+q, 47.32.Cc, 67.40.Vs, 67.57.Fg 



1. Introduction 

Recent experiments in 3 He-B 1 demonstrated that the fate of few vortices 
injected into a rapidly moving superfluid depends on dimensionless intrinsic 
temperature dependent parameter q, rather than on the flow velocity. At 
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q of order unity a rather sharp transition is observed between the laminar 
evolution of the injected vortices and the emerging turbulent many- vortex 
state of the whole superfluid. This added a new twist to the general theory 
of turbulence in superfluids developed by Vinen 2 ' 3 and others. Attempts to 
modify the theory in order to incorporate the new phenomenon, have been 
made in Refs. 5 and. 6 Here we follow the latter work 6 where we utilized the 
coarse-grained hydrodynamic equation for the dynamics of the superfluid 
with distributed vortices. In this equation the parameter q characterizes 
the friction force between the superfluid and normal components of the liq- 
uid, which is mediated by quantized vortices. According to this equation, 
turbulence develops only if the friction is relatively small compared to the 
inertial term, i.e. when q is below unity. We also argued that the developed 
turbulence must occur at q <C 1, and suggested that in this region, where 
the inertial term is dominating, there are at least two possible states of tur- 
bulence. One of them corresponds to the state discussed by Vinen, where 
the microscopic nature of quantized vortices is essential. While in the other 
state, turbulence does not depend on the circulation quantum, and thus the 
information on the underlying microscopic physics of quantum liquid is lost. 
Turbulence in this state becomes similar to turbulence in classical liquids, 
and it is also described by the Richardson cascade at least at the initial stage 
of the development. However, as distinct from the classical liquids, the final 
state of turbulence is determined not by viscosity but by the mutual friction 
parameter q. 

These differences from the classical turbulence arise because the hydro- 
dynamics of superfluid liquid exhibits new features with respect to conven- 
tional classical hydrodynamics, which become important when turbulence in 
superfluids is considered: 

(i) The superfluid liquid consists of mutually penetrating components 
- the viscous normal component and one or several frictionless superfluid 
components. That is why different types of turbulent motion are possible 
depending on whether the normal and superfluid components move together 
or separately. Here we are interested in the most simple case when the dy- 
namics of the normal component can be neglected. This occurs, for example, 
in the superfluid phases of 3 He where the normal component is so viscous 
that it is practically clamped to the container walls. The role of the normal 
component in this case is to provide the preferred heat-bath reference frame, 
where the normal component and thus the heat bath are at rest. The dissi- 
pation takes place when vortices move with respect to this reference frame. 
The turbulence in the superfluid component with the normal component at 
rest is here referred to as superfluid turbulence. We also assume that there is 
only a single superfluid component. For 3 He-B this means that we ignore the 
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spin degrees of freedom, assuming that all three superfluid spin components 
move together. 

(ii) The important feature of superfluid turbulence is that the vorticity 
of the superfluid component is quantized in terms of the elementary circu- 
lation quantum k (in 3 He-B k = irk/m, where m is the mass of 3 He atom). 
So, superfluid turbulence is the chaotic motion of well determined and well 
separated vortex filaments. 3 Using this as starting point we can simulate the 
main ingredients of classical turbulence - the chaotic dynamics of the vortex 
degrees of freedom of the liquid. However, to make the analogy useful one 
must choose the regime where the microscopic nature of the superfluid - the 
circulation quantum containing the Planck constant Ti and the mass of the 
atom m - is not involved. 

(iii) The most important distinction from the classical hydrodynamics 
is that the dissipation in vortex motion is not due to the viscosity term z^V 2 v 
in the Navier-Stokes equation. The superfluid component does not exhibit 
viscosity, instead the dissipation occurs due to the friction force acting on 
the superfluid vortex when it moves with respect to the heat bath (the nor- 
mal component). The force acting on a single vortex is proportional to the 
velocity of the vortex in the heat-bath frame. In the coarse-grained hydro- 
dynamics of the superfluid with the distributed vorticity, this leads to the 
force between the superfluid and normal components, which is proportional 
to their relative velocity v s — v n , and to the vortex density |V x v s |. As a 
result the dissipative term in the hydrodynamic equation is non-linear, and 
its structure reminds that of the non-linear inertial term. The relative mag- 
nitude of the two non-linear terms, the friction and inertial ones, is given by 
the dimensionless parameter q. Thus the quantity 1/q plays the role of the 
Reynolds number. This is an internal parameter of the liquid, as distinct 
form the Reynolds number in conventional liquids where it depends on the 
flow velocity and dimension of the system. When g<l, the inertial term 
is dominating, and this corresponds to a big Reynolds number in classical 
liquids. In this regime one can expect, that the turbulent state obeys the 
Richardson energy cascade, governed by the inertial term. We discuss here 
whether the Kolmogorov scaling law survives the non-canonical dissipation 
in superfhhds, i.e. whether or not the scaling is modified due to the non- 
linear dissipative term. 



2. Coarse-grained hydrodynamic equation 

The coarse-grained hydrodynamic equation is obtained from the Euler 
equation for the superfluid velocity v = v s after averaging over the vortex 
lines (see Refs. 7 and review 8 ). Instead of the Navier-Stokes equation with 
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the z/V 2 v term one has 

_ + V/i = v xu - (1) 

— a'(v — v n ) xw + awx(wx(v- v n )) . (2) 

Here v n is the velocity of the normal component; uj = V x v is the superfluid 
vorticity; uj = uj/lj; and the dimensionless parameters a' and a come from 
the reactive and dissipative forces acting on a vortex when it moves with 
respect to the normal component. These parameters are very similar to 
the Hall and conventional resistivities, p xy and p xx , in the Hall effect. For 
vortices in fermionic systems (superfluid 3 He and superconductors) they were 
calculated by Kopnin, 9 and measured in 3 He-B over a broad temperature 
range by Bevan et. al. 10 (see also, 11 where these parameters are discussed in 
terms of the chiral anomaly). We omitted the higher order terms, including 
the term discussed by Hall 7 which contains the bending energy of vortex lines 
(the condition for that is discussed in Section 6 where another dimensionless 
parameter will be introduced - the superfluid Reynolds number Re s = UR/k 
which contains the circulation quantum k). 

The terms in expression (1) are invariant with respect to the transfor- 
mation v — > v(r — ut) + u as in classical hydrodynamics. However, the 
terms in expression (2) are not invariant under this transformation, since 
there is the preferred reference frame in which the normal component is at 
rest. They are invariant under the full Galilean transformation when the 
normal component is also involved: v — > v(r — ut) + u and v n — > v n + u. 

The experiments in Ref. 1 were made in a rotating cryostat, where the 
normal component exhibits solid body rotation, v n = ft x r. However, here 
we study the local properties of turbulence and avoid this complication, 
assuming that v n is uniform. Actually this means that we consider the 
case of strong turbulence where the vortex density essentially exceeds the 
equilibrium vortex density in the rotating container: uj S> ft. We shall work 
in the frame comoving with the normal component, where v n = 0, but we 
must remember that this frame is unique. In this frame the equation for 
superfluid hydrodynamics is simplified: 

— + Vfi = (1 - a')v xw + awx(wxv). (3) 

After rescaling the time, t = (1 — a')t, one obtains an equation 

+ V/U = vxw + gwx(wxv), (4) 

which depends on a single parameter q: 
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Now in Eq. (4) the first three terms together are the same as the inertial terms 
in classical hydrodynamics. They satisfy the modified Galilean invariance: 

v(i, r) v(t, r - ut) + u . (6) 

In fact the transformation above changes the chemical potential, but this 
does not influence the vortex degrees of freedom which are important for 
the phenomenon of turbulence. In contrast, the dissipative last term with 
the factor q in Eq.(4) is not invariant under this transformation. This is in 
contrast to the conventional liquid where the whole Navier-Stokes equation, 
which contains the kinematic viscosity 

— +V/i = v x uj + »V 2 v , (7) 

is Galilean invariant, and where there is no preferred reference frame. 

This difference between the dissipative last terms in Eqs. (7) and (4) is 
very important: 

(1) The role of the Reynolds number, which characterizes the ratio of 
the inertial and dissipative terms in the hydrodynamic equations, is played 
by the intrinsic parameter 1/q in superfluid turbulence. This parameter does 
not depend on the characteristic velocity U and size R of the large-scale flow 
as distinct from the conventional Reynolds number Re = RU/u in classical 
viscous hydrodynamics. That is why the turbulent regime occurs only at 
1/q > 1 even if vortices are injected to the superfluid which moves with 
large velocity U. This rather unexpected result was obtained in experiments 
with superfluid 3 He-B. 1 

(2) In conventional turbulence the large-scale velocity U is always under- 
stood as the largest characteristic velocity difference in the inhomogeneous 
flow of the classical liquid. 12 In the two-fluid system the velocity U is the 
large-scale velocity of the superfluid component with respect to the normal 
component, and this velocity (the so-called counterflow velocity) can be com- 
pletely homogeneous (however, for the intermediate turbulent state obtained 
in the 3 He-B experiments in rotating vessel, the large-scale velocity field is 
inhomogeneous and we use U ~ QR, with Q being the rotation velocity and 
R - the radius of the container). 

(3) As distinct from the classical hydrodynamics, the energy dissipation 
is produced by the non-linear mutual friction (mf) term - the last term in 
Eq.(4): 

e mf = -E = -(v- ^X) = -q(\ • {Co x (Co x v))) ~ quv 2 . (8) 
ot 

This reflects the fact that the mutual friction depends on the velocity of the 
superfluid component with respect to the normal one, i.e. on the absolute 
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value of the velocity v in the frame of the container. Thus for a given value 
of vorticity u the largest dissipation occurs at the large-scale velocity U. On 
the other hand, according to the Kolmogorov scaling, the maximum vorticity 
uj is concentrated at the smallest scale. This gives the nontrivial dependence 
of the dissipation on scales: dissipation is not concentrated at the smallest 
scale only. 

(4) The onset of the superfluid turbulence was studied by Kopnin in 
Ref. 5 His Vinen-type model demonstrated that the initial avalanche-like 
multiplication of vortices leading to turbulence occurs when q drops below 
unity, which is in agreement with experiment. 1 The existence of two regimes 
in the initial development of vorticity is also supported by earlier simulations 
by Schwarz who noted that when a (or q) is reduced a crossover from a 
regime of isolated phase slips to a phase-slips cascade and then to the fully 
developed vortex turbulence occurs. 13 

(5) One can expect that the well developed turbulence occurs when 
q <C 1, and here we shall discuss this extreme limit. In 3 He-B the con- 
dition q <C 1 is realized at temperatures well below Q.QT C } However, we 
do not consider very low temperatures where instead of the mutual friction 
the other mechanisms of dissipation take place such as excitation of Kelvin 
waves 14 and vortex reconnection. 15 The latter leads to formation of cusps 
and kinks on the vortex filaments whose fast dynamics creates a burst of 
different types of excitations in quantum liquids: phonons, rotons, Kelvin 
waves and fermionic quasiparticles. A similar burst of gravitational waves 
from cusps and kinks on cosmic strings was theoretically investigated by 
the cosmological community (see e.g. 16 ), and the obtained results are very 
important for the superfluid turbulence at a very low temperature. 

(6) We expect that even at q <C 1, two different states of turbulence are 
possible, with the crossover (or transition) between them being determined 
by q and by another dimensionless parameter Re s = UR/k, where re is the 
circulation around the quantum vortex. The coarse-grained hydrodynamic 
equation (4) is in fact valid only in the limit Re s 3> 1, since the latter means 
that the characteristic circulation of the velocity r = UR of the large-scale 
flow substantially exceeds the circulation quantum re, and thus there are 
many vortices in the turbulent flow. When Re s is reduced the quantum 
nature of vortices becomes more pronounced. At some (still big) value of 
Re s we proceed from the type of superfluid turbulence which resembles the 
classical turbulence and is probably described by the (modified) Kolmogorov 
cascade, to the quantum regime which is probably described by the Vinen 
equations for the average vortex dynamics. 2 
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3. Kolmogorov cascade 



Let us first consider the possibility of a Kolmogorov- Richardson cascade 
in superfluid turbulence. Let us start with the analysis based on the dimen- 
sionalities and on the idea of the energy cascade. In the next Section we 
shall use a more detailed analysis based on model dynamical equations for 
the energy spectrum and on dispersion of dissipation over the scales. 

In classical turbulence, a large Reynolds number Re = UR/v S> 1 leads 
to well separated length scales or wave numbers. As a result the Kolmogorov- 
Richardson cascade takes place in which the energy flows from small wave 
numbers k m i n ~ 1/R (large rings of size R of the container) to high wave 
number ko = 1/ro where the dissipation occurs. In the same manner in 
our case of superfluid turbulence the necessary condition for a Kolmogorov 
cascade is a large ratio of the inertial and dissipative terms in Eq.(4), i.e. 



In the Kolmogorov-Richardson cascade, at arbitrary length scale r the 
energy transfer rate to the smaller scale, say r/2, is v^/t r , where t> 2 is the 
kinetic energy at this scale, and t r = r/v r is the characteristic time (turn- 
over time) . If there are no appreciable losses of energy in the intermediate 
scales, the energy transfer from one scale to the next must be the same for 
all scales. As a result one has 



This certainly occurs in classical liquids, and there is some indication that 
this might take place in superfluids too. 20 However, in superfluids the dissi- 
pation mechanism is different, which can change the overall picture: even if 
initially the energy flows to the smaller scale according to the Richardsonin, 
this does not mean that in the final steady state the Kolmogorov scaling will 
be established. 

In classical liquids, the energy dissipation due to viscosity IS 6 v isc ^ 
z^(Vjv) 2 . It is concentrated at the smallest scale ro, which is found from the 
energy balance e = e v i sc ~ uvq/tq = ve 2 / ? 'r 4 ^ 3 . This gives 



l/q > 1. 




(9) 



From this equation it follows that 



= e V3 r l/3 _ 



(10) 



r = l/fco ~ i?Re~ 3/4 . 



(11) 



In superfluid turbulence, the situation is more complicated. According 
to Eq.(8), in superfluids the dissipation is not a linear function of the velocity. 
In the non-linear theory, as follows from the Fourier analysis, contributions 
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from different scales are not independent. Since the vorticity uj r = v r /r is 
concentrated at the smallest scale, while the velocity - at the largest scales, 
one obtains that the overall dissipation due to mutual friction is e m f ~ qcooU 2 . 
Let us make a natural assumption that the Kolmogorov cascade stops when 
the overall dissipation becomes comparable to the energy transfer from scale 
to scale e in Eq.(9) 

e = e mf ~ q^U 2 ~ qU 2 ^ = qU 2 e^ r" 2 / 3 ~ ^ . (12) 

ro it 

Then, from the above equation one obtains the smallest length scale in the 
superfluid turbulence: 

r = 1/ko ~ q 3/2 R , (13) 

which depends on the internal parameter q instead of the Reynolds number 
Re in Eq.(ll) for classical liquids. The associated characteristic velocity vq 
and vorticity u)q at this scale are 

v ~ q l/2 U , uj = k v ~ ^ . (14) 

This consideration is valid when ro <C R and vq <C U, which means that 
1 is the condition for the Kolmogorov cascade. In classical liquids the 
corresponding condition for well developed turbulence is Re 3> 1. In both 
cases these conditions ensure that the kinetic terms in the hydrodynamic 
equations are much larger than the dissipative terms in a large enough iner- 
tial range. In the same manner as in classical liquids the condition for the 
stability of the turbulent flow is Re > 1, one may suggest that the condition 
for the stability of the discussed turbulent flow is 1/q > 1. This is supported 
by observations in 3 He-B where it was demonstrated that at high velocity U 
but at q > 1, turbulence is not developed even after vortices were introduced 
into the flow. 1 

In the Kolmogorov cascade both in classical liquids and superfluids the 
kinetic energy is concentrated at large scales comparable to the container 
size: 

E = I* -v 2 r = f R -(er) 2 / 3 = (ei?) 2 / 3 = U 2 . (15) 

Jro f Jtq f 

The dispersion of the turbulent energy in the momentum space is also 
the same as in classical liquid 

f R dr , rV R dk 9 fi/Rdke 2 / 3 „ 

E{k) = vl/k = e 2 / 3 fc- 5/3 . (16) 
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However, the dispersion of dissipation caused by the mutual friction is 
different from that in viscous classical liquids, where the main dissipation 
occurs at the smallest scale. This can change the whole pattern of the 
turbulent steady state. 

4. Dispersion of dissipation 

In superfluids the direct transfer of the kinetic energy to the normal 
component due to the mutual friction occurs at all scales. Let us consider 
the dispersion e™ f of this dissipation as a function of the length scale r. In 
the nonlinear term qojv 2 the given scale r contributes through the velocity 
field v and the vorticity field \oj\. The most important contributions are 
from the fluctuating vorticity u; r and the fluctuating velocity v r : 

ef = e « + e ( 2 ) = q(-)L0 2 . + q(uj)v 2 . (17) 

Id 

The term comes from the dissipation experienced by the fluctuating 
vorticity uj r at scale r. Assuming the Kolmogorov scaling one obtains the 
following estimation of this term 

4 1} = 9<-K 2 = q { ^ 2 . = ^ 2 . = e(r /rf 3 , (18) 

Here as before e = U 3 /R is the energy flux from scale to scale in the inertial 
range. We take into account that the vorticity moves with respect to the 
normal component with the characteristic velocity U and thus experiences 
the corresponding mutual friction. The Eq.(18) corresponds to the effective 
turbulent viscosity: = u e ffuj 2 with v e ff = qU 2 /ujq. As in the classical 
liquids this dissipation is peaked at the smallest scale where the vorticity is 
concentrated. 

The second term in Eq.(17) is the contribution of the velocity fluctu- 
ations to the friction experienced by the average vorticity (see Eq.(16) in 
Ref. 4 ). Under the same assumption of the Kolmogorov scaling, the second 
term has the following magnitude: 

4 2 > = qu v 2 = e(r/R) 2 / 3 . (19) 

Though the average vorticity is concentrated at the smallest scale, (lo) = ujq, 

(2) 

the term e r as a function of r is peaked at largest scale r = R where 
the velocity and kinetic energy are concentrated. This term which removes 
energy at large scales was also considered in classical liquids (see e.g. Eq.(5) 
in Ref. 17 with a = quo). 
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The other possible terms in the dissipation are smaller than these two 
if q <C 1. For example, one can add the term qv 2 co r . In the Kolmogorov 
cascade it does not depend on the scale r, and is always much smaller than 
the sum of the two main terms: qu r v 2 = qe <C q 2 / 3 e = min r (e^ + ef 2 ^). 
Thus at q <C 1 we have two peaks where the dissipation is concentrated. 
These two peaks at two extreme scales, tq and R, have equal magnitude if 
the Kolmogorov scaling is obeyed: 

e« ro ~ e^2 R ~ e = U 3 /R . (20) 

Let us now assume that dissipation which occurs at a given scale r is on 
the order of the losses of the kinetic energy at the same scale, i.e. dtE r ~ — e r . 
Since the dissipation is nonlinear, this is not the fact; nevetheless, this is a 
rather natural assumption. If one accepts this at least for the two scales, 
R and ro, then using the double-peak structure of the dispersion of the 
dissipation one comes to the following scenario of the cascade. The part of 
the initial kinetic energy of the superflow at large scale r = R is directly 
transferred to the normal component of the liquid due to the mutual friction 
in Eq.(19). The remaining comparable part of the flow energy experiences 
the Kolomogorov-Richardson cascade to the smaller scales until the next 
dissipation peak due to Eq.(18) is reached at r = tq. 

5. Possible modifications of the Kolmogorov cascade due to 
dispersion of dissipation 

If our assumption is correct, the above scenario with two peaks of dissi- 
pation does not change the Kolmogorov scaling discussed in Sec. 3, though 
some part of the kinetic energy is dissipated already at large scale. However, 
in principle, the dispersion of the dissipation can modify the energy cascade 
and the scaling law. Let us consider this using the following model of dissi- 
pation, which probably is not very realistic but allows us to find the modified 
scaling law within the model. Let us suppose that the two peaks overlap so 
that instead of the double-peak structure, the dissipation due to the direct 
transfer of the energy to the normal component is equally distributed over 
the scales: 

ef = e mf ~ qu U 2 , (21) 

and consider, whether the picture of the Kolmogorov-Richardson cascade 
survives or not under this new input. And if yes, how this modifies the 
Kolmogorov spectrum. Due to the direct losses of energies at each scale the 
energy transfer from scale to scale cascade d ecreases with decreasing scale, 
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so that the Eq.(9) does not hold any more: 

,,2 „3 rr3 
cascade u r u r / ~_ _ ('99 s l 

6r ~t r ~r R~ { > 

Instead we must write an equation for the r-dependent cascade - n ^ e s ^ eac jy 
state, which takes into account that the energy flux from the scale r to 
the neighbouring scale r/2 splits into the direct energy transfer from the 
scale r/2 to the normal component and the energy flux to the next scale: 
e cascadc _ £ cascadc _|_ Thus the difference between the energy fluxes at 
neighbouring scales is balanced by the direct energy losses at each step in 
Eq.(21) and one obtains the following equation for e c ascade : 

rd r e c r cadc = e mi , or rd r (^fj = e mf . (23) 

In terms of wave numbers k = 1/r and taking into account the time depen- 
dence, the energy balance in Eq.(23) can be written as 

dt{v 2 k ) + kd k (kv 3 k ) = -e mf , or 8 t E(k) + d k {k 5 / 2 E 3 / 2 ) = -^f . (24) 

where as before the energy E(k) = v\jk. 

We can also use the more sophisticated equations of Leith type de- 
scribing the diffusion of energy in the /c-space, which have been considered 
for turbulent cascades in conventional liquids. 18 For example, in a recent 
publication 19 the following diffusion equation was used: 

d t E = ^8 k {k ll ' 2 E l / 2 d k {E/k 2 )) + / - vk 2 E . (25) 

Here / is an external source of energy. In superfluids, instead of the last 
(viscous) term, there is a sink of energy caused by the direct transfer of 
energy to the normal component by mutual friction. In our model with 
constant dispersion of dissipation one has: 

dtE = ±0 k (k^E^d^E/k 2 )) + / - i £mf . (26) 

For this model the steady-state solutions of equations (24) and (26) are the 
same (assuming that the source term f(k) = at k > 1/R)- 

E(k) = C(e m{ f 3 k-VHn 2 / 3 ^ . (27) 

k 

They differ only by the prefactor C which is C = (24/1I) 2 / 3 for Eq.(26), and 
C = 1 for Eq.(24). 
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Eq.(27) demonstrates the main difference in the energy spectrum for the 
conventional and modified Kolmogorov cascades. In both cases the energy 
spectrum can be represented as 

E(k) = Ce 2 / 3 k- 5 / 3 F(k/k ) . (28) 

But in conventional turbulence one has F(k <C ko) = 1, while in the modified 
turbulence the function F(k/ko) logarithmically diverges at ko » k. 
From Eq.(27) it follows that the velocity 

v k = {kE)^~(e ml )^k-^\n^^ (29) 

monotonically decreases with k and approaches zero at k — ► ko. On the 
other hand the vorticity 

^, = ^~( emf )V3 fc 2/3 ln l/3| m 

first increases, then reaches its maximum 

Umax ~ (e m f) 1/3 fco /3 ( 31 ) 

at k max = ko/ \fe and finally almost abrubtly goes to zero within the same 
scale ko- 

The kinetic energy E{k) decreases first as /c~ 5 / 3 . Then, when the scale 
ko is reached, it vanishes at k = ko as (ko — fc) 2 / 3 . When k approaches 
ko, the non- linear inertial term in the hydrodynamic equation decreases and 
approaches the nondinear friction term. The latter becomes dominating and 
the cascade stops. 

Let us now find how the value of ko in Eq.(13) is modified by the log- 
arithmic function. The direct transfer of the energy to the heat bath e m f is 
determined by the maximal vorticity: 

e mf ~ quJmaxU 2 . (32) 

The total energy losses e is given by 

<- k ° dk 
'l/R 

As a result the energy spectrum becomes 



/ — e mf = e mf ln(fco-R) , (33) 
Ji/r k 



E(k) = Ce^k^^-f-. (34) 
In 1 ' 6 k R 
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From equations (31), (33) and (32) it follows that uj max ~ q^^k^XJ ~ 
e/{qU 2 ln(fco-R)). Then using the equation for energy at large scale R, U 2 = 
E(k = 1/R)/R, one obtains from Eq.(34) that the total energy losses to the 
heat bath is e = U 3 /R. Finally one comes to the logarithmic modification 
of equation (13) for the scale ko 

k = — , r = Rq 3/2 ln(k R) « Rq 3 / 2 ln{l/q) . (35) 

The characteristic velocity vo at this scale ro and vorticity uj max are: 

U u. 



7 ; 2 r . (36) 



/2gln(l/g) ' u fc 

Another model of the dispersion of the dissipation has been recently 
considered by Vinen in Ref. 4 Vinen took into account a single term in dissi- 
pation, Eq.(19), which is peaked at large scale. He solved the Eq.(26) with 
this type of dissipation and obtained the solution with a stronger modifica- 
tion of the Kolmogorov cascade. 

Here we discussed the steady state solutions of diffusion equation, but 
in principle, the time-dependent diffusion equation (25) for E(k, t) can be 
used for the analysis of the formation and decay of the turbulent state with 
different dispersion of the dissipation. 



6. Crossover to Vinen quantum turbulence 

Assuming that the superfluid turbulence is described by the non- 
modified Kolmogorov cascade, let us discuss the two possible regimes which 
occur at different ranges of q. At a very small q the microscopic nature of 
vortices with quantization of circulation becomes important. The condition 
of the validity of the coarse-grained hydrodynamic description used above 
is that the circulation relevant in the turbulent state can be considered as 
continuous. This means that the circulation at the smallest scale ro must be 
still larger than the circulation quantum, v^r® > k: 

^o r o = U Rq = q KKe s > k , Re s = . (37) 

K 

Thus the constraint for the validity of the Kolmogorov cascade is 

Re s > 4r > 1 • (38) 
q 2 

The same condition can be derived from the requirement that the character- 
istic scale ro must be much larger that the intervortex distance I. The latter 
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is obtained from the vortex density in the Kolmogorov state 

U 1 Re s 
nK = -~^ q = ■ (39) 

The condition I <C ro leads again to the equation v^r® > k and thus to the 
criterion (38). Actually the same condition allows us to neglect the term 
in the coarse grained equation one can neglect the term which takes into 
account the energy of the bending of vortex lines. 7 

In connection to Eq.(39) let us mention the experiments on 3 He-B in the 
rotating vessel. 1 There the turbulent state appears as the intermediate state 
between the initial vortex-free state and the final state where the superfluid 
experiences in average the solid-body rotation v = ft x r, where fi is the 
angular velocity of the rotating vessel. In the final equilibrium state the 
vortex density is ^-equilibrium = 2Q/k. From Eq.(39) and taking into account 
that £1 = U/R one obtains that in the intermediate turbulent state the 
vortex density must exceed the equilibrium density: 

«K ~ ^"-equilibrium • (40) 

The excess of vortices in the intermediate state has been observed in the 
experiment 1 and also in numerical simulations of this experiment. 21 

The criterion (38) contains the 'superfluid Reynolds number' Re s , which 
is determined by the microscopic characteristic of the vortex state - the 
circulation quantum k. This superfluid Reynolds number is responsible for 
the crossover or transition from the classical superfluid turbulence, where the 
quantized vortices are locally aligned (polarized) forming thick vortex tubes 
(Fig. 1 top), and thus the quantization is not important, to the quantum 
turbulence of isolated quantized vortex lines (Fig. 1 left) whose description 
was developed by Vinen. 

We can now consider the approach to the crossover from the quantum 
regime - the Vinen state which probably occurs when Re s q 2 < 1. According 
to Vinen 2 the turbulence in the quantum regime is characterized by a single 
scale. The characteristic velocity is the counterflow velocity U, while the 
characteristic length scale I is the distance between the vortices or the size of 
the characteristic vortex loops. It is determined by the circulation quantum 
and the counterflow velocity, / = Xk/U, where A is the dimensionless intrinsic 
parameter, which probably contains a' and a. The vortex density in the 
Vinen state is 22 

n v = r 2 ~ A 2 ^ = ^Re 2 , (41) 

It differs from the vortex density in the Kolmogorov state in Eq.(39) which 
depends not only on the counterflow velocity U, but also on the container 
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Fig. 1. Possible phase diagram of dynamical vortex states in the (Re s ,g) 
plane. At large flow velocity Re s 3> 1, the boundary between turbulent and 
'laminar' vortex flow approaches the vertical axis q = qo ~ 1 as suggested by 
experiment. 1 The thick line separates two regimes of superfluid turbulence 
occurring at small q: (i) the developed turbulence of the classical type, which 
is characterized by the Kolmogorov-type cascade possibly modified by the 
mutual friction; and (ii) the quantum turbulence of the Vinen type at even 
smaller q. 

size R. If the crossover between the classical and quantum regimes of the 
turbulent states occurs at Re s q 2 = 1, the two equations (41) and (39) must 
match each other in the crossover region. But this occurs only if A 2 ~ q. 
If A 2 ^ q there is a mismatch, and one may expect that either the two 
states are separated by the first-order phase transition, or there is an inter- 
mediate region where the superfluid turbulence is described by two different 
microscopic scales such as tq and I. 

In any case the crossover between the classical and quantum regimes 
occurs at the border of applicability of the coarse-grained equation (3), i.e. 
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we cannot use this equation to describe the Vinen state. For this state, in 
which the vortex lines are not locally aligned, the Vinen equations must be 
used. Note that there is a principle difference between these two turbulent 
states, which comes from the different dependence of the vortex density on 
U and R in Eq.(41) and Eq.(39). In particular, in the Vinen state, the 
mutual friction force between the normal and superfiuid components (the 
last term in Eq.(4)) obeys the cubic law suggested by Gorter and Mellink: 23 
/ ~ qnnyU ~ q\ 2 U 3 / k. In the Kolmogorov state, the friction force is 
obeying the quadratic law, / ~ qU ujq ~ U 2 /R. 

Based on the above consideration one may suggest the phase diagram 
of different regimes of collective vortex dynamics in Fig. 1. It is possible 
that the boundary between the 'laminar' and 'turbulent' regions reflects 
the process of developing of the turbulent vortex cluster: in the 'laminar' 
region turbulence cannot be started by a few injected vortices, while in the 
turbulent region the injected vortices lead to the vortex avalance and the 
turbulent state. 5 Actually this is what was observed in the experiment. 1 
However, the more fundamental role of this boundary, as the line of the 
phase transition between the vortex states, is not excluded: it is possible 
that in the laminar region the steady-state turbulence is simply unstable 
and decays. The instability of the turbulent vortex cluster in the lower part 
of the phase diagram was observed by Schwarz in his numerical simulations, 
but whether the boundary between the turbulent and laminar vortex flow 
approached the vertical axis q = qo ~ 1 at large velocity as suggested by 
experiment, 1 was not clear from the simulations. 22 



Superfiuid turbulence (as well as turbulence in classical 3-dimensional 
liquids) is a collective many-vortex phenomenon. Here we discussed the sim- 
ple case when the normal component is at rest, or its motion is fixed (which 
actually occurs in 3 He-B). Also we considered the dynamics of a single super- 
fluid component, i.e. we ignored the other possible superfiuid components 
in the multi-component superfluids, such as the spin degrees of freedom in 
3 He-B: spin currents and spin vortices. But even in this simple case there 
can be several different types of collective dynamical vortex states. Each of 
these vortex states can be characterized by its own correlation functions. For 
example, as the characteristic which distinguishes between different vortex 
states one can use the behavior of the loop function 



7. Discussion 




(42) 
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In the limit when the length L of the loop C is much larger than the inter- 
vortex distance I one may expect the general behavior g(L) ~ exp(— (L/Z) 7 ) 
where the exponent 7 is different for different vortex states. The asymptotic 
behavior of the loop function has been used for the description of the equi- 
librium phase transition in condensed matter and quantum field theory (see 
e.g. Ref. 24 > 25 ). It can be used also for the description of the non-equilibrium 
phase transition 11 where vortices also appear in the intermediate state ac- 
cording to the Kibble-Zurek mechanism. 26 ' 27 

In principle, the parameters a and a' may depend on the type of the 
dynamical state, since they are obtained by averaging of the forces acting on 
individual vortices. The renormalization of these parameters a{L) and ce'(L), 
when the length scale L is increasing, may also play an important role in the 
identification of the turbulent states, as in the case of the renormalization- 
group flow of similar parameters p xx and p xy in the quantum Hall effect (see 
e.g. 28 ). 

One can expect phase transitions between different states of collective 
vortex dynamics. One of such transitions, which appeared to be rather sharp, 
has been observed between the 'laminar' and 'turbulent' dynamics of vortices 
m superfluid 3 He-B. 1 It was found that this transition was regulated by the 
intrinsic velocity-independent dimensionless parameter q = a/(l—a'), rather 
than by flow velocity. However, it is not excluded that both dimensionless 
parameters a and a' are important. Moreover, it is also possible that only 
the initial stage of the formation of the turbulent state is governed by these 
parameters. 5 

Another transition (or maybe a crossover) is suggested here between the 
quantum and classical regimes of developed superfluid turbulence. We argue 
that there is a range of parameters in the classical region of the (Re s , q) plane 
where turbulence of vortex lines is described by the Kolmogorov-Richardson 
cascade, or by some modification of it caused by the dispersion of the dissi- 
pation. In this regime, turbulence is similar to that in classical liquids with 
modified dissipation. Thus the superfluid serves as a physically motivated 
example of the liquid with the non-canonical dissipation, which requires the 
general analysis of different models of dissipation and forcing (see e.g. 17 ). 

The turbulence in classical liquids is thought to be characterized by the 
dynamics of the vortex tubes, whose radii are of order of the dissipative 
Kolmogorov scale (see e.g. Ref. 29 and references therein). In superfluid tur- 
bulence, the radii of such tubes are determined by the parameter q acording 
to Eq.(35): vq ~ Rq 3 ^ 2 . In superfluids, the classical description is valid when 
these vortex tubes are fat enough, i.e. they contain many quantized vortices. 
The circulation around the tubes essentially exceeds the circulation quan- 
tum k, so that we can ignore the quantum nature of superfluid vortices. In 
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this region we can study analytically, numerically and experimentally all the 
phenomena which take place in classical liquids, including intermittency. 

The crossover to the quantum regime occurs when the circulation 
around the relevant vortex tubes becomes comparable to the circulation 
quantum k. According to Vinen, in the quantum regime, turbulence is 
represented by a single length scale - the intervortex distance I - and the 
single velocity scale U, which are related through the circulation quantum: 
I ~ k/U. It was suggested by Skrbek 30 that such a crossover between the 
quantum Vinen state and the classical Kolmogorov state of superfluid turbu- 
lence was probably observed in several experiments with the counterflowing 
superfluid 4 He. Skrbek identified these states with turbulent states I and II 
according to the Tough's classification scheme. 31 
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